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The ground state structural, electronic and magnetic properties of multiferroic TbMr^Os are 
investigated via first-principles calculations. We show that the ferroelectricity in TbMr^Os is driven 
by the non-centrosymmetric magnetic ordering, without invoking the spin-orbit coupling and non- 
collinear spins. The intrinsic electric polarization in this compound is calculated to be 1187 nC- 
cm -2 , an order of magnitude larger than previously thought. 
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Multiferroics with magnetic and electric ordering 
united in a single phase were thought to be rare 

a a. 

Surprisingly, a large class of manganese oxides (RMn03 
0, |||, and RMn 2 5 0, & with R= Y , Tb , D y, etc ) 
has recently been discovered to be multiferroic. Unlike 
the traditional multiferroics where the two order param- 
eters only couple weakly [l[, the newly discovered mate- 
rials possess strong magnetoelectric (ME) coupling, re- 
sulting in various novel physical effects. One of the most 
prominent examples is given by TbMn205[5|, [8[, which 
displays clear correlated anomalies of the dielectric con- 
stant e with the spin ordering [5]. More strikingly, the 
electric polarization in this material can be reversed by 
applying a magnetic field [5[ . The remarkable ME effects 
revealed in these materials have attracted great attention 
EL 0, H 0, 0, 0, 0, 0, Ed, El El because of the fascinating 
physics and their potential applications in novel multi- 
functional ME devices. 

Although great effort has been devoted to understand- 
ing the fundamental mechanism of the giant ME cou- 
pling, our knowledge of the manganese oxide multifer- 
roics is still very limited and full of puzzles. For example, 
experimental data show that the structure of TbM^Os 
has space group Pbam [13], which includes spatial inver- 
sion (R -1 ) symmetry. It is therefore puzzling that the 
compound develops spontaneous electric polarizations. 
It has been suspected |6|, [9[ that the actual symmetry 
group of TbMn 2 05 is Pb2\m, allowing polarizations. In- 
deed, there are several experiments supporting this hy- 
pothesis Nevertheless, no direct evidence of 
the lower symmetry has yet been found jf| 0]. Theo- 
retically, the origin of the giant ME coupling and the 
ferroelectricity in these materials is still under intensive 
debates |( 



|, [ill, [12[ . One of the fundamental questions 
that remain unsolved is whether the spin-orbit interac- 
tion []j], EH is essential for the ME coupling and ferro- 
electricity in these materials. 
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In this letter, we carry out a first-principles study of 
TbMn205, to identify its crystal structure and clarify 
the microscopic origin of the ferroelectricity. We com- 
pute physical quantities that can be directly compared 
to the experiments. To the best of our knowledge, no 
such study has yet been done for TbM^Os (and other 
RMn205 compounds [7]), because it has very compli- 
cate incommensurate anti-ferromagnetic (AFM) struc- 
ture with the propagation vector k « (0.48,0,0.32). To 
accommodate the magnetic structure, one needs a huge 
supercell, which is computationally prohibitive. Instead, 
we use a 2x1x1 supercell, equivalent to approximating 
the propagation vector k = (0.5, 0, 0). The validity of this 
approximation will be justified later in the text. Our re- 
sults show that the ferroelectricity in TbM^Os is driven 
by the non-centrosymmetric magnetic ordering and the 
asymmetric exchange interactions, without invoking the 
spin-orbit coupling and non-collinear spins. The intrin- 
sic electric polarization in this compound is calculated to 
be 1187 nC- cm -2 , much larger than previously thought 
for this type of compound [l|, [5[ . 

The calculations are based on the density functional 
theory (DFT) within the spin-polarized generalized gra- 
dient approximation (GGA) [14J implemented in the Vi- 
enna ab initio Simulations Package (VASP) [15|,ll6[. The 
projector augmented- wave (PAW) pseudopotentials [13] 
with a 500 eV plane-wave cutoff are used. A 1 x 2 x 4 
Monkhorst-Pack k-points mesh converges very well the 
results. We use the collinear spin approximation without 
the spin-orbit coupling. Our results agree very well with 
the known experiments, indicating that these approxima- 
tions capture the essential physics in TbM^Os. 

The crystal structure of TbM^Os is orthorhom- 
bic, with four TbM^Os formula units (32 atoms in 
total) per primitive cell, containing Mn 4+ 06 octahe- 
dra and Mn 3+ Os pyramids [13]. We relax the crys- 
tal structure beginning with the experimental structural 
parameters [Hj, listed in Table U The most stable struc- 
ture we found has the spin configuration identical to that 
was proposed in Ref. |g, as illustrated in Fig. [TJa). It also 
has an energetically degenerate structure, shown in Fig. 
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TABLE I: Comparison of the calculated and measured struc- 
tural parameters of TbM^Os. The lattice constants are 
given in A. Atoms that occupy the same Wyckoff positions 
are shown only once. 



Theory (Pb2im) Experiment (Pbam) 

a b c a b c 

a, b e 7.3014 8.5393 5.6056 7.3251 8.5168 5.6750 

Tb a+ 0.1410 0.1733 0.1399 0.1726 

0.6404 0.3270 

Mn 4+ 0.0001 0.5003 0.2558 0.5 0.2618 

Mn 3+ 0.4012 0.3558 0.5 0.4120 0.3510 0.5 
0.9016 0.1456 0.5 

01 0.0008 0.0002 0.2709 0.2710 

2 0.1645 0.4480 0.1617 0.4463 
0.6648 0.0517 

3 0.1560 0.4329 0.5 0.1528 0.4324 0.5 
0.6571 0.0655 0.5 

4 0.3977 0.2077 0.2438 0.3973 0.2062 0.2483 
0.8959 0.2919 0.7579 



(a) (b) 




FIG. 1: (Color online) The ground state spin configurations 
for two energetically degenerate structure L and R. The 
diamonds and triangles denote the Mn 4+ 06 octahedra and 
Mn 3+ 05 pyramids respectively. The dashed, single and dou- 
ble lines represent J3, J a and J5 exchange interactions respec- 
tively, following Ref. @. 

[2(b). We denote the two structures "left" (L) and "right" 
(R) respectively. In these magnetic structures, Mn 4+ 
form an AFM square lattice in the ab plane, whereas 
Mn 3+ couples to Mn 4+ either antiferromagnetically via 
J4 along a axis or with alternating sign via J3 along b 
axis. Mn 3+ ions in two connected pyramids also couple 
antiferromagnetically through J5. Here, we adopt the 
notations J3, J4 and J5 from Ref. [y, and define the J3 to 
be the Mn 4+ - Mn 3+ superexchange interaction through 
pyramidal base corners, and J4 the superexchange inter- 
action through the pyramidal apex, as indicated in Fig. 
[TJ We label the two different Mn 4+ chains along the a 
axis I, II respectively, also following Ref. 6. The mag- 
netic structure of R can be obtained from L by shifting 
chain II to the right (or to the left) by one unit cell along 
the a axis [6[. 

The calculated structural parameters for the structure 
L are listed in Table U comparing with the experimen- 
tal data, whereas the structure R is a mirror image of 
L about ac-plane. The calculated structural parameters 
are in extremely good agreement with the experimen- 



TABLE II: Comparison of the Mn 3+ positions in Pbam sym- 
metry and Pb2±m symmetry. 





Pbam 


Pb2 1 m 


h 


(x, y, 1/2) b(l) 


( x, y, 1/2) 


h 


(-x+1/2, y+1/2, 1/2) b(l) 


(-x+1/2, y+1/2, 1/2) 


h 


(x+1/2, -y+1/2, 1/2) b(2) 


(x+l/2+fe, -y+l/2-iy, 1/2) 


h 


(-x, -y, 1/2) b(2) 


(-xrtx, -y-Sy, 1/2) 



tal data. The errors of lattice constants are about 1%, 
typical errors for GGA. The atom positions are also ex- 
tremely close to what was obtained experimentally. How- 
ever, the small atomic displacements lower the structural 
symmetry to the long searched Pb2\m polar group. To 
see how this happen, we take Mn 3+ ions as an example. 
In the Pbam structure, Mn 3+ has one Wyckoff position 
(h in Wyckoff notation) that has four equivalent sites 
shown in Table Hfl However, in the Pb2\m structure, it 
splits into two Wyckoff positions b(l), b(2), each has two 
equivalent sites. Sx and Sy in Table HT1 denote the atomic 
displacements from the high symmetry positions along 
the a and b axes respectively. The displacements along 
the a axis are of mirror symmetry, whereas the displace- 
ments along the b axis are not, allowing polarizations. 
The atomic displacements can be easily calculated from 
Table H We see that the displacements are extremely 
small, usually are of the order of ~ 10 -4 of the lattice 
constants and the largest atom displacements come from 
Mn 3+ , Sy ~ 10 -3 of the lattice constants. Therefore the 
low symmetry structure can not be directly determined 
experimentally, and only the anomalies of the atomic dis- 
placement parameters (ADPs) were observed [6]. We 
also artifically construct a high symmetry structure by 
symmetrizing structure L and R according to the Pbam 
symmetry, which we refer as structure H in the following 
discussions. 

The calculated density of states (DOS) of structure 
L (and also R) is shown in Fig[2l The DOS for spin- 
up and spin-down electrons is identical as expected for 
an AFM state. The DOS exhibits a small but clear band 
gap (about 0.4 eV), confirming the experimental fact that 
TbMn205 is an insulator. However, it is well known that 
GGA greatly underestimates the band gap, especially for 
the 3d compounds. The local magnetic moments are es- 
timated for Mn 3+ to be ~ 2.37 and for Mn 4+ to be 
~ 1.64 [1b, in good agreement with the refined magnetic 
moments [6[. 

One of the strong evidences suggesting that the low 
temperature structure of TbM^Os has space group 
Pb2\m is that some Raman active modes become also 
IR active [10| at low temperature, which is forbidden by 
higher Pbam symmetry. To clarify this problem, we an- 
alyze the zone-center phonons. The symmetry analyses 
are performed on the 32-atom primitive cell [18]. For the 
high symmetry structure (Pbam), the total 96 modes, are 
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FIG. 2: The total DOS for TbMn 2 5 of structures L and R. 
The dashed line indicates valence-band maximum. 



decomposed into 8 irreducible representations (irreps): 

T = B lu © B 2u © B 3u © A3 © B lg ®B 2g ®B 3g ®A u . (1) 

Among them £? 3n , B 2u and B lu modes are IR active, po- 
larized along the a, b and c axes, respectively. B\ g , B 2gi 
B 3g and A g modes are Raman active while A u modes are 
silent. As we see, the Raman- and IR- active modes do 
not couple. However, the real crystal structure possesses 
a Pb2\m symmetry, whose 96 phonons can be decom- 
posed into 4 irreps: 



T = A 1 © B 1 © B 2 © A 2 . 



(2) 



We found all modes are Raman active, among them Ai, 
Bi and B 2 are also IR active, and are polarized along the 
6, c, and a axes respectively [10j . Detailed analyses show 
that the A\ modes are coupled from the A g and B 2u 
modes. We then calculate the phonon frequencies via 
a frozen-phonon technique 11811 . as well as the oscillator 
strengths for the IR modes [19[. The calculated phonon 
frequencies are all in good agreement with the experi- 
ments, generally within 20 cm -1 from the experimental 
values. All phonons are found to be stable, i.e., no soft 
phonon has been found as in traditional ferroelectrics. In 
this letter, we show only the b- axis-polar A\ phonons in 
Table Hill with frequencies and their oscillator strengths. 
To make a good contact with experiments, we divide the 
phonons into two presentations A g and B 2u according 
to their major symmetry character. The experimental 
frequencies of the A g phonons are extracted from Ra- 
man spectra [20], while the B 2u phonon frequencies are 
extracted from IR spectra [10j. The overall calculated 
oscillator strengths of A\ modes are in good agreement 
with the experimental values [lOj] . More importantly, the 
693 cm -1 A g Raman mode, that was found also IR ac- 
tive with oscillator strength 5\ =0.001 in the experiment 
is well reproduced in the calculations, with oscil- 
lator strength 5^=0.004, therefore, confirming that the 
ground state structure is indeed of Pb2\m symmetry. 



TABLE III: Calculated phonon frequencies (oo) and oscillator 
strengths (Sx) of IR- active Ai modes compared with experi- 
mental values. The modes are divided into E>2u and A^ repre- 
sentations according to their major symmetry character. The 
experimental values of E>2 W modes are taken from Ref. HoL 
whereas those of A g modes are taken from Ref. I20I . 
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Experimentally observed IR active and 703 cm 1 measured by 
Aguilar et al.flOj 



We next calculate the spontaneous polarization in the 
compound using the Berry-phase theory of polarization 
implemented in VASP [21]. The intrinsic polarization in 
this material is calculated to be 1187 nC-cm -2 along the 
b axis. This value is an order of magnitude smaller than 
that of the traditional ferroelectrics. e.g., BaTiOs, yet 
it is about 30 times larger than the currently measured 
experimental value (~ 40 nC-cm -2 [5]) for this com- 
pound. The large discrepancy between calculated and 
experimental polarizations might come from the approx- 
imations we used in the calculations. For example, we 
approximate the magnetic propagation vector k z =0.32 
by zero. We also ignore the spin-orbit coupling and as- 
suming collinear spins. Without these approximations, 
the polarization might be smaller. On the other hand, 
the experiment [5[ measured polycrystalline samples, in 
which grains polarize along different directions, cancel- 
ing each other, therefore might greatly underestimated 
the intrinsic polarization. We believe a high quality sin- 
gle crystal sample should enhance the measured electric 
polarization. 

To further elucidate the origin of the polarization, we 
also calculate the spontaneous polarization for the high 
symmetry structure H and get P= 228 nC-cm -2 . In 
this case, the crystal structure has Rr 1 symmetry, how- 
ever, because the special spin configuration (see Fig. [T|) 
of the Mn 3+ -Mn 4+ -Mn 3+ chains along the b axis breaks 
the R~ x and the combined spatial inversion and time re- 
versal [(RT) -1 ) symmetry of the magnetic structure, the 
electron wavefunctions have lower symmetry than the lat- 
tice, resulting in polarization. The electronic symmetry 
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FIG. 3: (Color online) An asymmetric- spin- chain model 
showing the sketch of energy surfaces of structure L and R vs. 
the atomic displacements from the high symmetry structure. 
The diamonds and triangles denote the Mn 4+ 06 octahedra 
and Mn 3+ 05 pyramids respectively. 

breaking will further couple to the lattice and lead to lat- 
tice distortion. Therefore, when holding atoms fixed at 
the centrosymmetric structure, turning on the magnetic 
order does two things: it generates a purely electronic 
polarization of 228 nC-cm -2 , and it also applies forces 
to the atoms. These forces give rise to atomic displace- 
ments that yield an additional 959 nC-cm -2 . 

We now discuss the microscopic mechanism of the cou- 
pling between magnetic order and lattice. Figure [3] de- 
picts the energy surfaces of the structure L and R vs. 
atomic displacements. The two structures degenerate 
in energy at the high symmetry structure H (Pbam). 
The magnetic structures of L and R can be simplified 
as Mn 3+ -Mn 4+ -Mn 3+ segments linked along the b axis 
via J5. Inside the segments, Mn 3+ interact with Mn 4+ 
via exchange interaction J3. In each segment, the Mn 3+ 
ions above and below Mn 4+ ions have opposite spins. For 
structure L, Mn 4+ always has the same sign of spin with 
the lower Mn 3+ , but opposite sign with the upper one in 
each segment, whereas the opposite is true for structure 
R. In the high symmetry structure H, the Mn 3+ and 
Mn 4+ ions with the same sign of spin could move closer 
to each other to minimize the exchange energies [8]. The 
formal description of the coupling between the spin and 
the lattice is given by [H, [l9[ 

ij\ 

where Jij is the exchange interaction between the mag- 
netic moments S^, Sj of the z-th and j-th atoms, and ua 
is the A-th zone-center phonon calculated for structure 
H. This term is not zero, provides the magnetic struc- 



ture (spin configuration) {S^} does not have the R~ x and 
the {RT)~ X symmetry, as e.g., in TbM^Os. Therefore 
the high symmetry structure H is an unstable point on 
the energy surface as illustrated in Fig. [3l The structure 
will spontaneously relax to L or R according to the spin 
configuration at the high symmetry point. Furthermore, 
as long as the asymmetric spin structures of the Mn 3+ - 
Mn 4+ -Mn 3+ chains are persevered, the propagation vec- 
tor k z will not change the essential physics in this system 
[22l | which justify our approximation of setting k z =0. 

To conclude, we have shown via a first-principle study 
on TbMn205, how the ferroelectricity is driven by the 
non-centrosymmetric magnetic ordering, without invok- 
ing the spin-orbit coupling and non-collinear spins. We 
believe that this work sheds new light on the funda- 
mental mechanism of the giant magnetoelectric coupling 
in the multiferroics, especially for the manganese ox- 
ides. Surprisingly, the calculated intrinsic polarization 
in TbMn205 is as large as 1187 nC-cm -2 , much larger 
than previously thought for this type of compound. The 
mechanism revealed here also provides useful guidance 
for searching for novel magnetoelectric materials. 
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